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General  Abstract: 

Nicrostructures  were  studied  in  the  presence  of  generalized  defects, 
both  point  repulsive  potentials  and  slit  openings  in  barriers.  A  number  of 
transport  effects  were  discovered.  For  isolated  defects  a  resonance  was 
found  to  occur  at  the  threshold  energies  of  a  quantum  well.  A  scaling 
relation  was  found  for  scatterers  near  single  barriers.  A  rich  spectrum  of 
effects  was  found  for  single  defects  near  and  within  double-barrier  struc¬ 
tures.  The  electronic  analogue  of  single  slit  diffraction  was  shown  to 
exhibit  high  fringe  visibilities. 

Written  presentations  of  research  results: 

"Electronic  Transport  in  Nicrostructures  with  Defects,"  (Master's 
Thesis  of  Brent  Haukness)  -  This  presents  a  detailed  description  of  the 
methods  used,  with  62  figures  (143  pp.).  It  is  available  on  microfilm  from 
University  Microfilms  Inc.,  Ann  Arbor,  Michigan.  Abstract: 

The  effects  of  elastic  point  scatterers  on  the  electronic  transport 
properties  of  microstructures  are  examined  using  a  generalized  transfer 
matrix  approach.  Using  a  delta-function  potential  to  model  a  defect, 
transmission  characteristics  are  calculated  for  a  number  of  different 
microstructures.  Among  these  is  the  double  barrier  resonant  tunneling 
diode.  A  number  of  effects  are  noticed.  These  effects  depend  sensi¬ 
tively  on  the  position  of  the  defect  relative  to  specific  features  of 
the  potential  structure,  and  a  saturation  of  the  effects  is  noticed 
with  increasing  defect  strength.  Some  of  the  peculiarities  of  the 
delta-function  defect,  are  also  examined,  and  a  comparison  is  made  to 
defects  modeled  with  a  Gaussian  potential. 

"Analysis  of  Electron  Diffraction  in  a  Novel  Field-Effect  Transistor," 
by  A.  N.  Kriman,  G.  H.  Bernstein,  B.  S.  Haukness  and  D.  K.  Ferry.  -  This 
original  paper  was  published  in  Superlattices  and  Microstructures,  Vol.  6, 
No.  4,  pp.  381-366,  1969,  and  presented  as  a  poster  at  the  1988  Conference 
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on  Superlattices,  Microstructures  and  Microdevices  held  in  Trieste,  Italy. 
It  provides  a  theoretical  analysis  of  the  behavior  of  the  QUADFET  device. 
(The  vork  of  G.  H.  Bernstein  and  D.  K.  Ferry,  and  the  work  of  A.  M.  Kriman 
related  to  this  project,  were  not  paid  by  this  contract.  The  vork  of 
Haukness  was.)  This  paper  represents  calculations  related  to  the  main  topic 
of  geometric  effects  in  scattering  in  a  number  of  ways:  (i)  It  provided  an 
"entry  level"  problem  for  the  student  (Haukness)  which  taught  him  some  of 
the  skills  that  he  was  to  need  in  mastering  the  main  project,  such  as  using 
the  CONVEX  computer  with  UNIX  operating  system  and  with  local  plotting 
facilities.  (ii)  It  involved  a  microstructure  in  two  dimensions  with  non¬ 
trivial  geometry.  The  situation  treated  -  a  long  thin  barrier  with  a  slit, 
is  in  fact  a  nontrivial  microstructure  with  a  defect.  The  "defect"  is  a 
slit  which  intentionally  built  into  the  barrier.  (iii)  The  mathematical 
approach  was  closely  related  to  that  used  to  solve  the  general  problem 
ultimately  treated  by  Haukness.  Abstract: 

Ve  analyze  a  field  effect  transistor  whose  operation  utilizes  the 
quantum  diffraction  of  electrons  by  a  narrow  slit  in  the  gate. 
Calculations  are  performed  using  a  low-energy  conformal  mapping  tech¬ 
nique  which  does  not  assume  small-angle  scattering.  It  is  shown  that 
the  device  will  exhibit  far-field  quantum  diffraction  similar  in  ap¬ 
pearance  to  the  Fraunhofer  patterns  observed  in  optics.  The 
diffraction  pattern  can  be  detected  as  current  collected  at  a  number  of 
narrow  Schottky  contacts  which  together  comprise  a  "viewing  screen". 
Fringe  visibilities  on  the  order  of  0.5  are  predicted.  A  number  of 
applications  of  the  device  are  discussed. 


"Geometric  Effects  of  Scattering  in  Microstructures,"  by  A.  M.  Kriman, 

R.  P.  Joshi,  B.  S.  Haukness  and  D.  K.  Ferry.  -  This  original  paper  will  be 
published  in  Solid  State  Electronics,  Vol.  32,  No.  12,  pp.  1597-1601,  1989,  „r 
and  was  presented  as  a  poster  at  the  1989  Conference  on  Hot  Carriers  in 
Semiconductors  held  in  Scottsdale,  Arizona.  Abstract: 

Trans  Eer  Matrix  techniques  are  used  to  study  elastic  scattering  by  a. 
point  defects  embedded  in  quasl-one-dimensional  microstructures.  This 
makes  possible  an  exact  analysis  of  phenomena  that  arise  from  breaking  T 
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of  transverse  translation  invariance.  The  dependence  of  transmission 
probability  on  scatterer  position  is  studied  for  parallel  transport  in 
quantum  veils  and  for  perpendicular  transport  across  single  and  double 
barrier  structures.  It  is  found  that  in  laterally  confined  structures, 
delta>function  and  other  extremely  sharp  models  of  a  single  defect  lead 
to  sharp  resonances  when  such  defects  are  veil  isolated.  Such  features 
are  associated  vith  multiple  reflections  betveen  the  lateral  confining 
potential  and  the  defect  potential.  In  single-barrier  structures  vith 
a  single  nearby  defect,  an  approximate  scaling  behavior  is  found  that 
relates  transmission  for  defects  at  different  distances  to  that  at  a 
fixed  distance  vith  different  energy  scales.  In  double  barrier  reson¬ 
ant  tunneling  diodes  (DBRTDs),  the  position  of  the  transmission  peak  is 
affected  primarily  by  defects  vithin  the  quantum  veil  region.  The 
height  of  the  transmission  peak  is  very  sensitive  to  the  positions  of 
defects  vithin  that  region,  acting  essentially  as  a  probe  of  the 
resonance  vave  function.  Defects  in  front  of  a  DBRTD  also  affect  the 
valley  current  by  modifying  the  longitudinal  component  of  the  incident 
momentum. 

Educational  accomplishments: 

One  graduate  student,  Brent  S.  Baukness  (a  U.S.  citizen),  vas  hired  for 
the  year  1988.  As  a  result  primarily  of  the  vork  that  he  performed  under 
this  contract,  he  earned  a  Haster's  degree  in  Electrical  Engineering.  Re 
had  offers  to  attend  graduate  school  tovard  his  Ph.D. ,  but  he  preferred  to 
enter  industry;  he  accepted  an  offer  from,  and  has  been  vorking  since 
January  for  Texas  Instruments,  Inc. 

It  vas  not  found  possible,  as  envisaged  in  the  original  budget,  to  hire 
another  graduate  student  vith  adequate  qualifications  to  perform  vork  re¬ 
quired  under  the  contract.  Instead,  Ravindra  P.  Joshi,  vho  had  recently 
received  his  Ph.D.  in  Electrical  Engineering  from  Arizona  State  University, 
vas  hired  on  a  part-time  basis.  (He  vorked  full-time  over  all,  as  he  also 
vorked  part-time  for  R.  0.  Grondin.)  This  vork  constituted  post-doctoral 
training  alloved  under  his  visa,  and  it  did  significantly  Increase  the 
breadth  of  his  experience  beyond  vhat  he  had  learned  as  a  graduate  student. 
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EFFECTIVE  POTENTIAL  FOR  MOMENT-METHOD  SIMULATION  OF  QUANTUM  DEVICES 
A.  M.  Kriman,  J.-R.  Zhou,  N.  C.  Klufcidahi,  H.  H.  Choi,  and  D.  K.  Ferry 


Center  for  Solid  Stale  Electronics  Research 
Arizona  Slate  Universiiy,Teinpe,  Arizona  85287-6206 


ABSTRACT 

In  the  simulation  of  submicron  devices,  complete  quantum  descriptions  can  be  extremely  compuutionally  intensive,  and 
reduced  descriptions  arc  desirable.  One  such  description  itfilizes  a  few  low-order  moments  ol  the  momentum  distribution  that 
are  defined  by  the  Wigner  function.  Tvra  major  mfftculiies  occur  in  applying  this  moment  method:  (i)  An  independent 
calculation  is  required  to  find  quantum  mechanically  accurate  initial  conditions,  (ii)  For  a  system  in  a  mixed  state,  the 
hierarchy  of  lime  evoluiion  equations  for  the  moments  does  not  close.  We  diescribe  an  approach  to  solve  these  problems. 
The  initial  distribution  is  determined  in  equilibrium  by  means  of  a  new  effective  potential,  chosen  for  its  ability  to  treat  the 
sharp  potential  features  which  occur  in  helerostructures.  It  accurately  describes  barrier  penetration  and  repulsion,  u  well  as 
quantum  broadening  of  the  momentum  distribution.  The  moment  eouaiion  hierarchy  is  closed  at  the  level  of  the  second- 
moment  lime  evolution  equation,  using  a  closure  that  is  exact  for  a  shifted  Feimi  distribution.  Band-bending  is  included  by 
siinulianeous  self-consistent  deierminatian  of  all  the  moments. 


KEYWORDS 

Submicron  device  simulation;  quantum  barrier  repulsion  and  penetration;  moment  method. 


INTRODUCTION 

The  development  of  nanometer-scale  electronic  devices  has  required  the  modeling  of  quantum  mechanical  effects  in  device 
simulations.  A  number  of  such  simulations  have  been  develof^  for  the  resonant  tunneling  diode  (RTD),  whose  operation 
depends  essentially  on  the  purely  quantum  mechanical  effect  of  barrier  tunneling.  Because  the  Rlio  is  a  quasi-one- 
dimensional  device,  its  analysis  is  numerically  tracublc,  even  when  "self-consistency"  (electrostatic  fields  due  to  mobile 
charge  redistribution)  is  included  (Kluksdahl  ei  of.,  1989).  More  general  devices,  whose  active  regions  cannot  be  treated  in  a 
oiK-dimensianal  approximatkin,  present  a  computatianal  challenge. 

These  computational  proUeitu  already  occur  in  the  simulation  of  classical  (non-quantum  mechanical)  devices.  Qassical 
^vices  are  adequately  described  by  the  single-paiiicle  distribution  function  which  gives  the  density  in  the  six-dimensional 
phase  space  of  momentum  and  space  coordinates.  Near  equilibrium,  the  momentum  distribution  at  each  point  is  well 
approximated  by  a  shifted  eouilibrium  distribution.  This  observation  leads  to  the  successful  and  widely  used  "moment 
method"  (BHHekJn,  1966  and  1970).  In  this  approach,  instead  of  determining  the  full  distribution  function  as  a  function  of 
six  arguments,  one  uses  a  reduced  description  based  on  a  few  momentum  integrals  -  moments  of  the  distribution  function  - 
as  functkms  of  spatial  poridon.  The  moments  obey  time  evolution  equations  ("imment  equations")  that  are  integrals  of  the 
Boltzmann  equation. 

The  quantum  mechanical  generalization  of  die  classical  distribution  function  is  the  Wigner  distribution  function,  a  partial 
Fourier  transform  of  the  dr^y  matrix,  lafrate,  Grubin  and  Ferry  (1981)  developed  a  quantum  generalization  of  the  moment 
method,  based  on  integrals  of  the  Wigner  function.  Interestingly,  the  first  three  quanhim  moment  equations  are  formally 
identicaiiotheclassiaJones.  Thus,  ifouannmiefrectt  are  to  be  Mluded  in  the  simulaiion,  then  they  must  already  be  present 
in  the  initial  conditioos.  This  fact  it  also  known  from  simulationt  based  directly  on  the  Wigner  function  (Camthen  and 
Zachariaaon,  1983):  For  motion  fat  a  parabolic  potential,  the  full  classical  and  quantum  dittribiition  functions  obe^  identical 
equations  of  motion,  so  the  diflerenoe  between  cuosied  and  quantum  behavion  mutt  be  "contoined”  in  the  initial  distribuiion. 
While  this  demonttrMet  that  tome  aspectt  of  the  quaniuitVcIasticai  distinction  mutt  be  contained  in  the  initial  conditkmt,  it  is 
nevetthelen  the  case  that  for  motioo  HI  general  potentiids,  the  classicsl  and  quantum  evohnkin  equations  also  differ. 

The  moment  equations  form  a  hierarchy,  with  the  time  derivative  of  each  moment  determined  in  part  by  the  gradient  of  the 
next  moment  far  the  hientchv.  In  a  practical  calculation,  one  computes  only  a  few  Cow-order  moments,  closing  the  hierarchy 
by  making  tome  approximaaon  for  the  highest  moment  used.  Sucnanappmachit  justified  by  the  fact  that  at  equilibrium,  the 
exact  ditmbution  it  completely  detctminM  by  the  density  and  temperature  (the  zeroth  and  second  momentt).  Away  from 
equilibrium,  the  clcnure  approximation  corresponds  to  an  auumption  regarding  the  form  of  a  quasi-equilibrium  distribution, 
la  our  simulations,  we  use  the  first  three  momentt  (number,  momentum,  and  energy  densltiet).  Our  doaute  it  initially 
derived  by  auuming  that  third  momems  are  smaU.  but  it  it  exactly  true  for  a  large  dus  of  disadbuiion  functions,  induding 
the  thiAed  Feimi-Dwae  dittifliution. 

In  one-dimcational  RTD  models,  it  it  possible  to  compute  a  compleiefe  seif-conslttent  Wigner  distribution  function  in 
equilibrium  by  performing  a  theiriwlly  wdghtcd  turn  over  all  states.  This  hat  provided  an  iniiid  distribution  which  it 
accurately  stationary  uadcr  numeried  integratioa  of  the  Wigner  funedoo  time  evolutiaa  ei  radon.  In  contrast,  a  three- 
dimenaiond  device  has  many  mom  sinist,  and  if  He  SehrMHngcr  equation  it  not  tepambi a  is  difficali  even  to  find 
individad  states,  let  alone  pmorm  a  cotnpleta  sum.  Another  method,  based  on  path  Integra.'^  does  nrovlda  a  tystamatic 
approach,  but  it  wo  is  conipiitadonally  intensive,  and  has  so  date  been  applied  ptiinarily  to  one -dimensional  prebiw.  |See 
Mason  and  Hess,  (1989),  Regitier  at  d.  (I9M).  and  refemneas  therein.)  What  It  thoefora  requited  it  a  cowputatlotially 
efficiem  way  of  determining  an  equilibrium  (iaidal)  distribution,  or  set  of  momems,  thm  would  cotnplemeni  the 


oominiiukirally  eflicienl  moment  method  At  noted,  the  moments  to  detennined  mutt  mwiifett  quantum  mechanical  effecu: ' 
quantum  mechanical  conectiont  cannot  ante  independently  from  the  time  evolution  equations. 

Some  of  the  authon  (Kiiman,  Zhou  and  Feny.  1989)  have  lecently  proposed  a  tutittical  (“eflective  potential”)  method.  This 
provides  a  ouantum  mechanically  accurate  equilibrium  distribution  which  we  will  use  to  sadtf^  the  leouirement  for  accurate 
mitial  conditions.  The  effective  potential  dekribes  hairier  penetration  and  repulsion,  along  with  ^uiiibrium  broadenin|  of 
the  momentum  disiribiition  which  does  not  occur  classically.  Other  effective  potentials,  doived  rigorously  u  semklassical 
approximations  of  the  paidtion  function,  cannot  be  implied  in  the  presence  to  typical  device  models.  Inde^,  these  previous 
efiwtive  potentials  diverge  everywhere  if  the  potentiu  is  infinite  in  any  finite  region.  The  new  effective  potential  is  a  well- 
behaved  nmctional  of  the  teal  potendal,  and  hu  the  proper  asymptotic  Mhavior  near  an  infinite  potential  burier,  giving  enors 
in  density  at  the  few  per  cent  level.  Like  previous  effective  potentisls,  it  also  approaches  the  coirect  high-tempnature  limits, 
in  the  foUowing,  we  will  recapitulate  the  method  and  outline  its  derivation. 


EFFECTIVE  POTENTIAL 

We  find  an  initial  condition  in  equilibrium  by  using  a  recently  developed  effective  potential  method  (Kriman,  Zhou  and  Ferry, 
1989).  This  method  is  similar  in  form  to  earlier  statistical  approximations  develop  by  Giachetti  and  Tomtti  (i98S;  1986) 
and  ^  Kleineit  (1986)  with  Rynman  (1986).  In  common  with  these,  the  new  method  computes  the  density  classically  from 
an  effective  potential  which  incorporates  a  local  average  of  the  actual  potential.  The  local  average  rroresents  the  quantum 
position  uncertainty  of  a  particle  whose  energy  -  and  thus  momentum  -  is  not  completely  uncertain.  The  information  about 
the  energy  distribution  anses  fiom  the  finite  tempaature,  and  thus  the  range  over  which  the  position  must  be  averaged  is  of 

the  order  of  the  thermal  wavelength  Xth"  A  V|V2m  . 

The  main  respects  in  which  earlier  effective  potential  methods  differ  from  the  one  we  use  here  have  to  do  with  the  effect  of 
infinite  potentials.  The  Giachetti-Tognetli  arid  Feynman-Kleinett  effective  potential  approaches  usume  that  the  underlying 
potentisJ  is  smooth,  and  proceed  from  there  to  derive  analytically  approximations  which  ate  exact  in  the  high-temperature 
limit.  This  is  supported  te  the  fact  that  in  the  high  teinperaiure  limit,  the  lowest-order  quantum  correction  to  the  density  is  of 
relative  order  P(a.th)^  v  ,  where  V*  is  the  second  oerivative  of  the  potential  (Husimi,  1940).  Many  of  the  potentials 
commonly  used  in  modeling  solid  slate  -  and  particularly  hclerostniciure  -  devices,  contain  potentials  that  are  not  smooth. 
The  extreme  case  of  this  is  the  infiniie  wall  poiuiiial,  used  where  appropriate  to  set  the  density  to  xero. 

The  previous  effective  potentials  all  diverge  for  this  potential,  u  they  do  for  any  potential  whose  integral  in  any  finite  region 
is  innniie.  This  is  because  the  effective  potentials  are  linear  averages  of  the  ttu  potential,  and  the  weighting  in  space  falls 
smoothly  to  sero.  The  new  potential  method  that  we  have  introducixl  uses  nonlinear  averages,  and  remains  well-defined  in 
the  presence  of  infiniie  barrios. 

The  nonlinear  effective  potential  is  derived  in  the  spirit  of  density  functional  theory  (DFT)  (Jones  and  Gunnarsson,  1989). 
The  DFT  is  based  on  the  faa  (Hohenberg  and  Kohn,  1964:  Levy,  1979;  Mermin,  IMS)  that  the  election  energy  is  a  unique 
functional  of  the  electron  density.  Although  the  exact  functional  is  unknown  and  presumably  quite  co^lex  (Lieb,  1982),  it 
is  possible  to  develop  approximate  functimials  that  are  adjusted  to  be  quite  accurate  for  a  set  eSf  simplified  cases  (Kohn  and 
Sham,  1965).  Anak^ously,  the  nonlinear  effective  potential  is  based  on  the  observation  that  the  density  is  (by  constniction)  a 
unique  functional  of  the  potential.  While  the  exact  formal  expressian  of  this  functional  can  be  difficult  to  evaluate,  it  is  again 
possible  to  develop  approximate  functionals  by  adjusting  pariunem  in  a  suitable  Ansatz.  We  have  also  used  certain  simple, 
exact  properties  of  die  density  to  constrain  the  form  of  the  effective  potential.  These  constraints  were  that  (a)  the  effective 
potenaal  approach  the  real  one  at  high  temperature,  and  that  the  effective  potential  preserve  (b)  honnogeneity  and  (c) 
separability  propetiies  of  the  real  potential. 

The  high-iemperatnte  oonstraim  stems  from  one  motivation  of  effective  potentials  in  general:  that  in  the  classical  regime  the 
density  is  given  simply  bv  an  exponmtial  of  the  teal  potential,  and  that  this  regime  is  approached  at  high  temperatures.  One 
simple  coose(|uenoe  6t  this  constraint  is  that  various  weighting  functions  must  be  normalised.  Another  is  that  the  local 
averaging  must  occur  over  length  scales  that  decrease  with  increasing  temperature.  By  dimensional  analysis,  for  a  sharp 
potential  feature  this  length  scale  must  be  essentially  the  thermal  wavelength. 

The  homogeneity  constraint  we  impose  is  that  a  constant  change  far  the  real  potential  lead  to  an  equal  change  in  the  effective 
potentiaL  This  is  a  suprisingly  shingent  constraint  For  exam^,  consider  die  potential  has  the  general  functional  form 

V.(*)»h;'(<h.(V)>(x)),  (1) 


-t 


srbere  ha  is  a  function  with  inverse  lb  under  compositian: 


f(Jt)nh;'(h.(fXJr)).  (2) 

Here  h  is  a  real  function  of  a  real  argument,  and  nonlocality  It  introduced  by  the  average  It  can  be  proven  that  only 
two  fbtnis  of  h  ate  allowable  in  (I)  under  the  constraint  of  homogeneity.  Tmm  are  hi.(v)  ■  v,  which  leads  to  the  ordinary 
aridunctic  (Uncar)  averago,  and  h^v)  ■  cxp(-)^),  which  leads  to  the  exponential  averages  diM  we  mutt  ultimately  use  if 
htfinitepoientiabaioioMadmisaaole.  (Slight^  mote  general  k'tae  allowable,  but  lead  to  the  tame  averages.) 

The  avarabtlily  oonurahH  we  use  it  that  if  a  potential  it  snaraUe  in  Canetian  coordinatet,  then  its  effective  potential  win  be 
Mi^ns^is^i  s^i^l^tr^i^il^t,  ^^T^tls  prooerty  Mlowtfmn  the  factottsability  of  the  density  when  the  pMMdal  it  separable.  Theweak 
condition  Imposed  it  aoi  tufnclent  to  imply  diat  aeparability  will  be  conserved  for  the  other  coorriinaie  tyttemt  in  which  the 
SchrOdiager  equation  it  teparable.  However,  by  meant  of  a  standard  argument  from  the  kinetic  thco^  of  gates,  it  does 
follow  that  avnaget  neh  u  the  one  appearing  n  (I)  most  be  defined  by  a  weighting  function  that  is  uauttbn.  (See,  for 
tuafk,  (CatltBan,  1971)). 
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Numerical  paramelen  in  the  effective  potential  weie  set  by  optiimaingihe  fit  to  the  density  in  semi-infinite  wall  potential, 
whose  density  is  known  exactly  (Kriman,  Kluksdahl,  Feiry,  1987).  llie  temperatuK  and  mass  can  be  scaled  out  of  this 
potential,  so  the  results  have  a  very  general  validity.  After  requiring  that  the  feeing  asymptotic  behavior  of  the  density  near 
the  wall  be  given  exactly  (which  in  itself  rales  out  many  functional  Kmu  for  the  effective  potential),  the  maximum  fractional 
error  was  imnimized.  (jioosing  the  form 


•O 

fJw(Y/p)V^x), 

(3) 

with 

w(s) »  Wt  exp(-tu7>8s) , 

(4) 

we  found  that  wiy«l.4427,  a<«14.9442,  and  5-3.82dS.  The  weight  in  the  average  defining  is  a  Gaussian  of  rms  width 
2o»3.3220kTH(‘l^)*^-  It  should  be  understood  that  these  four  parameters  ate  the  result  of  a  two-parameter  optimization:  (i) 
The  small-x  asymptotic  behavior  of  the  density  (p~x^)  requires  w(0)«(in2)~'.  (ii)  The  high-temperaiur^classical-limit 
constraint  requites  w  to  be  normalized,  so  S  depencB  on  a  in  a  way  that  has  nothing  to  do  with  the  model  potential  chosen  for 
optimization.  The  peak  in  the  distribution  w  occurs  when  o»0.733X'ni.  Thus,  the  effective  potential  represents  an  average  of 
the  real  potential  over  a  range  comparable  to  the  Ibennal  wavelength.  (This  was  expected,  but  it  was  not  imposed  as  an 
independent  condition.) 

We  used  this  effective  potential  to  compute  the  density  in  equilibrium  of  a  resonant  tunneling  diode.  The  diode  consists  of 
0.3cV  barriers,  Snm  long,  separated  by  a  quantum  well  region  Snm  long.  The  election  is  assumed  to  have  the  a  mass  0.067 
times  the  free  eleettan  mass,  appropriate  for  GaAs.  In  Fig.  I  the  result  trf  the  effective  potential  calculation  is  compared  with 
the  exact  numerical  result  (Kluksdw,  et  al.,  1989).  comnuted  by  summing  over  scattering  eigenstates. 


Fig.  1.  Electron  density  from  exact  approach  (sum  over  scattering  states,  full  line) 
and  effective  potential  (broken  line). 


MOMENT  MEmOO 

A  reduced  description  of  semiconductor  systems  that  has  been  successful  for  classical  problems  is  the  moment  method 
(BKhekjcr,  1966  and  1970).  in  ouanium  proMems,  the  analogue  of  the  distribution  functioiK  which  provides  a  "complete" 
single-particle  description  is  the  wigner  aistribution  function  (Brittin  and  Oiappell,  (1962);  Camiihers  and  Zachanason, 
(1983);  Balazs  and  Jennings,  (1984);  Hillery,  <l  rtf.  (1984)).  It  is  a  Fourier  transform  of  the  density  matrix  p(x,z')  ■ 
(tyl(z)t|r(z‘)>  (tyl  and  tpare  TteM  operaioa],  taken  with  re^tect  to  the  difference  coordinate; 

•• 

|rfy  p(x+ Jji-f)  .  (5) 

The  Wigner  distribation  fiinctioti  fw  is  the  Weyl  transform  of  the  classical  single  panicle  phase  space  density.  That  is,  it  is  a 
particular  nuantum  generaliaiion  of  the  phase  qrace  density.  It  is  a  real  function,  and  one  can  immediately  define  local 
tnomems  of  the  momentum  distribution  by 
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(6) 


»!,(*)■  |<fpfwUiP.>)P" 

related  to  local  averages  of  the  momeiuum  and  its  powers  by 

0>">-l*i/P.  (7) 

where  p  ■  Mo-  It  is  possible  to  derive,  from  the  exact  quantum  mechanical  ^uation  of  motion  of  the  Wigner  function,  a  set  of 
moment  equations  analogous  to  the  moment  equations  of  classical  semiconductor  problems  (lafrate,  Gnibin  and  Ferry, 
1981).  The  first  three  moments  evolve  in  lime  acrarding  to 


where  the  inhomogenems  terms  on  the  right  contain  the  eflects  of  collisions.  One  may  reasonably  expect  that  for  a  good 
approximation  of  a  sMlisiical  system's  behavior.  It  would  suffice  to  evolve  just  the  first  three  elemenu,  corresponding  to 
dmsiiy,  momentum,  and  energy.  However,  as  noted  earlier,  closing  the  hierarchy  at  an^  finite  order  requires  some  extra 
information,  typically  in  the  form  of  an  expression,  using  lower-order  moments,  of  the  highesl-oider  moment  appearing  in 
the  e^uatioiis.  In  perticular,  the  last  of  our  nroment  equations  involves  the  gradient  of  113,  for  which  rw  lime  evolution 
equation  ia  given.  One  wav  of  "closing"  the  moment  equations  is  to  make  some  approximation  within  the  derivation  of  the 
moment  equations,  in  this  implicit  approach,  a  sepanie  closure  expression  may  not  occur,  and  the  approximation  is  built  into 
the  time  evolution  equaiioru,  which  are  also  usually  referred  to  as  ("the")  moment  equations.  To  remove  any  ambiguity,  let 
us  sute  that  where  we  refer,  in  this  paper,  to  the  formal  equivalence  between  the  first  three  classical  and  quantum  moment 
equations,  we  are  referring  spedfinlly  to  (8),  the  moment  ^uations  derived  without  any  hierarchy-closing  Ansaiz. 

In  the  special  case  of  a  pure  state,  the  "statistical"  system  is  described  completely  by  two  teal  quantities;  f  >  Ile(y)  -t- 
flttsfy).  Thus  one  expects  that  it  should  be  possible  to  close  the  imment  equatioru  at  the  level  of  die  two  moments.  In  fact, 
lafrate  and  coworkeis  (1981)  have  shown  that  for  a  pure  sute, 

•*1  *7  a* , 

The  second  term  is  (up  to  a  multiplicative  factor)  the  "quantum  potential"  of  Bohm,  and  using  this  closure  in  the  first  two 
equations  of  the  hierarchy  (8)  yields  the  same  moriKnis  as  a  direct  solution  of  the  Schrbdinger  equation. 

For  the  general  (mixed  sute)  case,  we  first  derive  our  closure  expression  by  regarding  the  skewness  of  the  momentum 
distribution  as  sniall.  We  write 

and  find  that  (fiT)  m  (p2)  .  (p)2,  mg  of  course  (8) «  0.  If  we  make  the  assumption  of  approximately  Oaunian  distributions, 
or  some  other  weaker  assunnption  which  causes  the  ritewness  {fi)  to  vanish,  then  we  can  conclude  that 

(pJ)-3<p*)<p>-2(p)J.  (9) 

This  closure  rule  has  a  very  useful  property:  It  is  easily  demonstrated  that  a  large  class  of  distribudons,  including  the 
common  drifted  Maxwellian  and  drifted  Fermi  distributions,  satisfy  this  relation.  In  general,  it  bolds  for  any  distribudon  of 
the  form  f(x^,f)  -  I^,(p-(p(*,<)»7,i). 

It  is  impoitani  to  nme  that  while  the  form  of  our  closuic  is  ((HP»^  -  0  in  the  classical  case,  different  closuies  arise  when 
sitnilar-appeating  approximations  are  made  in  different  contexts.  In  panicular,  Flosujc^  and  Rboa^-Biown  (1985)  and 
Siroscio  (1986)  use  an  br^icii  ctosure  that  is  essendally  equivalem  to  setting  -  ((Mp))  v  ■■  0  and  solving  the 

moment  equation  for  113  (the  lowest-oirier  equation  not  shown  in  (8)1.  The  closuie  obtained  in  this  wav  does  contain  a 
quamuffi  correction.  How^,dwteu  in  most  pcttatbativcaptaaches,  increasing  orders  of  h  are  assodaied  with  increasing 
orders  of  derivatives  of  dw  potcntiaf,  making  it  inaiyropriaie  for  the  problems  we  wish  10  neat 


SUMMARY 

We  have  proposed  m  approach  10  die  momem  mediod  for  quantum  devices,  using  an  elTecdve  poienda)4iased  density  which 
is  desigiM  to  mimic  tavortam  feaiaies  of  the  me  density,  such  u  tunneling,  barrier  penetration,  and  the  complementary 
barrier  repulsion.  The  choice  of  poremial  functional  is  guided  by  the  rcquiremem  to  sadsfy  a  number  of  fundamenul 
consirainis.  We  describe  an  explicil  closure  for  the  moment  cquadons  at  the  level  of  the  kinetic  energy  time  evolution 
equation  which  is  sadsfled  exactly  ty  ail  "drifted"  distributions. 
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We  analyze  a  field  effect  transistor  whose  operation  utilizes  the  quantum  diffraction  of  electrons  by  a  narrow  slit  in 
the  gate.  Calculations  are  performed  using  a  low-energy  conformal  mapping  technique  which  does  not  assume 
small-angle  scattering.  It  is  shown  that  the  device  will  exhibit  far-field  quantum  diffraction  similar  in  appearance  to 
the  Fraunhofer  patterns  observed  in  optics.  The  diffraction  pattern  can  be  detected  as  current  collected  at  a  number 
of  narrow  Schottky  contacts  which  together  comprise  a  "viewing  screen”.  Fringe  visibilities  on  the  order  of  0.5  are 
predicted.  A  number  of  applications  of  the  device  are  discussed. 


1.  Introduction 

Refinements  in  nanolithography  and  epitaxial  growth 
techniques  have  made  possible  devices  with  progressively 
smaller  features,  allowing  the  observation  of  retd-space  wave 
properties  of  electrons  in  the  solid  state.  Paradigmatic 
examples  of  these  are  resonant  tunneling  diodes  [I]  and 
small  ring  structures  [2]  exhibiting  Aharonov-Bohm  effects 
[3].  The  signature  of  the  wave  mechanical  effects  is 
typically  some  non-monotonic  feature  in  a  plot  of  current 
versus  electric  or  magnetic  field.  Such  features  arise  as  the 
interference  between  different  electron  paths  alternates 
between  predominantly  constructive  and  predominantly 
destructive. 

In  at  least  two  respects,  the  devices  that  have  been  made 
or  suggested  thus  far  generally  test  only  special  cases  of 
interference  phenomena.  First,  with  only  one  exception  that 
we  are  aware  of  (Furuya  (4)),  there  have  been  no  diffraction 
devices,  such  as  gratings  or  slits,  in  which  "bright"  and 
"dark”  regions  may  be  observed  simultaneously.  While 
electron  diffraction  has  been  demonstrated  and  utilized  [S] 
for  measurement,  the  diffraction  pattern  is  generally  detected 
a  macroscopic  distance  from  the  diffracting  structure,  rather 
than  within  the  microscopic  region  that  produces  it. 

Second,  the  interference  is  usually  between  quantum 
paths  which  differ  in  some  discrete  way,  as  in  an  Aharonov- 
Bohm  ring,  rather  than  among  a  continuous  range  of  paths, 
as  in  sin^e-slil  diffraction.  In  a  path  integral  view,  for 
instance,  the  resonant  tunneling  diode  operates  by  the 
interference  between  paths  differing  in  the  number  of  times 
they  traverse  the  quantum  well.  In  the  Aharonov-Bohm 
rings,  the  interfering  paths  differ  discretely  in  their  winding 


number  around  the  flux  enclosed  by  the  ring.  In  the  device 
proposed  by  Furuya,  diffraction  is  effected  by  a  transverse 
potential  grating  in  the  base,  and  the  paths  are  distinguished 
by  the  opening  of  the  base  through  which  they  cross. 

We  describe  a  device  which  differs  fiom  existing  devices 
in  both  of  the  above-described  respects.  This  "QUADFET" 
(Quantum  Qiffraction  EED  [6]  is  essentially  a  single-slit 
diffraction  experiment  realized  within  the  two-dimensional 
elect^  gas  of  an  inverted  HEMT  layer.  The  analog  of  the 
viewing  screen  is  a  set  of  forward-biased  Schottky  contacts, 
or  "collectors",  which  serve  as  the  drain  (see  Fig.  ]).  The 
slit  width  can  be  varied  by  adjusting  the  gate  potential,  and 
lobes  of  the  diffraction  pattern  are  detected  as  variations  in 
the  current  reaching  different  collectors. 

It  is  important  to  consider  the  survival  of  the  diffraction 
pattern  in  the  presence  of  collisions.  Beenakker  et  al.  [7] 
performed  elegant  experiments  which  show  that  elastic 
collisions  can  be  ignored  over  path  lengths  of  at  least  3  pm. 
Their  results  demonstrated  not  only  that  the  phase  memory 
of  the  electrons  was  conserved,  but  also  that  the  elastic 
collisions  did  not  cause  the  paths  to  deviate  significantly.  If 
elastic  collisions  did  cause  significant  deviation  of  the  path, 
then  the  Fraunhofer  diffraction  pattern  in  the  QUADFET 
would  be  washed  out  by  the  collisions  within  an  elastic  mean 
free  rath  length. 

This  paper  is  concerned  primarily  with  how  the 
QUADFfc  l  would  work  in  the  absence  of  scattering.  A 
more  central  problem  is  whether  the  device  can  be  made  to 
work  given  the  fundamental  differences  between  the  usual 
optical  difffaction  experiment  and  the  electronic  device  we 
have  proposed. 
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Fig.  1 .  Top  view  of  the  QUADFET. 


The  first  difference  is  in  the  source:  it  is  difficult  to 
monochromate  and  focus  electrons  in  a  semiconductor 
device,  whereas  c^tical  devices  benefit  from  readily  available 
laser  sources  shining  through  transparent,  nondispersive 
media.  While  some  monochromation  is  possible  by  the 
interposition  of  a  gate,  in  the  manner  of  the  THETA  device 
[8],  this  brings  into  play  the  constraint  of  statistics.  Because 
electrons  are  femions,  their  density  in  phase  space  is 
bounded  (by  2/h3  on  a  coarse  grain).  Thus,  focussing  and 
monochromation  diminish  the  current,  so  it  is  useful  to 
design  a  device  which  functions  well  with  an  unfiltered 
source  distribution. 

It  is  clear,  therefore,  that  we  must  deal  with  a  spatially 
extended  source  having  a  distribution  of  energies.  These 
give  rise  to  space-  and  time-domain  "coherence"  effects 
which  have  been  studied  in  ordinary  light  optics  (9,10J.  In 
the  specific  examples  that  are  usually  considered,  however, 
some  kind  of  small-angle  approximation  is  usu^ly  made. 
Without  some  such  approximation,  the  usual  Huygens- 
Fresnel-Kirchoff  approach  becomes  complicated. 
Therefore,  in  analyzing  the  QUADFET  we  have  developed  a 
new  approach  related  to  those  used  in  waveguide  analysis. 

2.  Device  Model  and  Diffraction  Amplitudes 

Our  centra]  problem  is  to  give  an  adequate  treatment  of 
the  quantum  states  in  the  vicinity  of  the  slit.  Initially, 
therefore,  we  solve  the  Schrbdinger  equation  for  a  potential 
that  consists  of  a  single  slit  gate  lying  along  the  y  axis: 

V('x,y)  =  VoCf)  5(x),  (1) 

(x  is  the  coordinate  perpendicular  to  the  gate,  within  the 
two-dimensional  electron  gas),  where  we  consider  the 
special  case 


0  forlyl<<V2 
oo  forlylia/2 


(2) 


The  gate  is  assumed,  like  the  corresponding  optical  barrier, 
to  be  infinitesimally  thin  (i.e.,  «  8(x)).  This  is  justified  if 
the  Fermi  wavelength  is  long  compared  to  the  actual  gate 
length.  The  barrier  "strength"  Vo(y)  is  nonnegative,  so  there 
are  no  bound  states.  Thus,  a  basis  of  eigensutes  for  this 
potential  consists  of  scanering  states  incident  from  the  right 
and  left.  Completely  general  momentum-normalized  [11] 
left-incident  states  can  be  written  (using  continuity  of  the 
wavefunction  at ;(  =  0)  in  the  form  [  1 2] 


‘l'k(x,y)  =  -  sin(i:;jX)  exp(ii:j,y)  e(-x)  +  Ea  .  (3) 


where  0(tt)  =  {  1  if  u  >  0,  0  if  u  ^  0  )  is  the  unit  step 
function,  and 


’^ix.y)  = 


2jt 


<y  F(qyJcy,E)  exp(/qxU:l+i<7>y) .  (4) 


(The  integration,  here  and  wherever  explicit  bounds  are  not 
given,  is  over  the  whole  real  line.)  This  eigenstate  describes 
an  electron  incident  with  wave  vector  k  =  (t„  ky)  (where  t* 

>  0).  The  sine  term  in  (3)  is  the  standing  wave  that  arises  if 
there  is  perfect  specular  reflection  atx  =  0.  The  difference 
between  this  and  the  exact  wavefunction  is  Za,  which 
contains  components  scattered  into  final  wave  vectors  q  = 
(9i>  Qy)-  Because  our  chosen  potential  (2)  reflects 
specularly  almost  everywhere,  the  correction  Za  S®**  to 
zero  quicUy  with  increasing  distance  from  the  slit,  and  F  is  a 
smooth  function  of  q,.  For  a  permeable  (V()J[y)  <  «<>)  or 
rough  gate  potential,  F  would  have  a  component 
proportional  to  &(qy-ky),  implying  finite  probabilities  for 
transmission  and  nonspecular  reflection. 

The  electron  generally  is  scattered  with  transverse  wave 
vectors  taking  the  full  range  of  values  -«>  <  <  »•.  The 

tot^  energy  E  =S  ^k^l2m  is  conserved  (i  =  Ikl),  so  q^  =  k\ 
which  m^es  q^  an  implicit  double-valued  function  of  qy  for 
a  given  energy.  The  sign  of  q*  must  be  chosen  to  satisfy  the 
scattering  boundary  conditions:  £|,  should  have  waves 
propagating  away  from  the  gate  (q^  >  0  when  lqj,l  <  *),  or 
exponentially  damped  away  from  the  slit  (-iq,  >  0  when  lq,l 

>  k).  Note  that  although  they  do  not  contribute  to  the 
current,  these  exponenti^ly  damped  waves  are  an  essential 
part  of  the  state,  and  contribute  to  the  density  near  the  gate. 
(By  reflection  symmetry,  the  right-incident  states  have  an 
expresision  similar  to  Eq.  (3).) 

The  scattering  state  (3)  is  completely  specified  by  F.  An 
earlier  study  [12]  concerned  with  idealized  contacts  treated  a 
geometrically  similar  problem.  Generalizing  the  method 
developed  there  (see  in  particular  the  case  treated  in  the 
Appendix),  one  can  write  down  an  implicit  solution  for  a 
potential  of  the  general  form  (1):  F  is  that  function  which 
satisfies 


dq,  F(q, ,*,;£)  (q,  +  ^o(y))  exp[i(qj,-*j,)y]  =  1  (5) 


for  all  y. 

In  the  earlier  study,  it  was  found  that  the  dominant 
behavior  of  transmission  amplitudes  (to  q.  from  ky)  arose 
from  the  factor  in  (4),  while  the  energy-dependence  of  F 
was  a  higher-order  effect.  In  order  to  keep  the  problem 
tractable,  therefore,  we  have  assumed 


F(qy,ky;E)  s  F(q,,*,;0)  s  F(qy,ky) .  (6) 

This  approximation  is  completely  analogous  to  one  made 
in  calculations  of  microwave  waveguide  characteristics  [13] 
and  of  fluid  flow  around  obstacles  and  through  constrictions 
[14].  The  computational  advantage  of  this  approach  is  the 
possibility  of  applying  conformal  mapping  techniques,  as  we 
shall  do.  We  apply  a  mapping  into  elliptic  coordinates 
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H  s  w  =  cosh-'(2iz/a) ,  (7) 

with  z  =  X  +  iy. 

The  infinite  bairiers  of  the  gate  now  coincide  with  ^=0  (jt 
=  0,  y  <  -a/2)  and  %  =  it  {x  =  0,  y  >  ail).  We  can  now 
implement  the  potential  simply  as  a  Neumann  (4*  =  0) 
boundary  condition,  and  we  can  expand  in  the  complete, 
discrete  basis  of  separated  states  that  results: 

Ik  =  X  sin(B^)  e'^,  (8) 

f|sl 


where  the  superscript  labels  the  half-plane  of  x:  c,(±)  =  c„W 
for  0  <  ^  <  71  (x  >  0),  and  c,(-)  for  7i  <  §  <  2ji  (x  <  0). 
These  expansion  coefficients  are  implicitly  functions  of  k  = 
[iky,  ky);  and  can  be  determined  by  continuity  conditions  at 
the  sht  (the  locus  x  =  0,  ly  I  <  a/2  is  just  n  =  0).  Continuity 
of  the  wavefunction  immediately  implies  =  C,<->  s  C„. 
Denvative  continuity  takes  the  form 

-|j-H'(p=0+,4)  =  +|j.>i'(n=0-^,27t-^) . 

Inserting  (3),  we  find 

c„  =  -r' Ukyafl) ,  (9) 

K 


where  is  the  n^-order  Bessel  function  of  the  first  kind 
[15],  By  identifying  (4)  and  (8),  and  taking  Fourier 
components,  we  find  that 

FidyJcy)  =  ni-))”  Jnikya/l)  UQyafZ) .  (10) 

We  will  be  concerned  with  diffraction  patterns  in  the 
Fraunhofer  regime:  current  will  be  collected  many 
wavelengths  away  from  the  slit.  As  noted  earlier,  F  is  a 
smooth  function  of  Qy,  so  we  can  use  a  stationary  phase 
analysis  of  (4)  in  this  regime  to  find  that  the  wavefunction  is 
asymptotically 

'I'i,(x=rcos0,y=rsin0)  ~ 

(L  'y  — 

- — ^  F(ksin0,A:y)  COS0  exp[r(kr- j)]  . 

in  polar  coordinates;  the  expression  is  valid  at  large 
distances.  More  precisely,  it  is  valid  whenever  'Ikir  cos0 
is  much  smaller  than  the  Qy  -scale  for  variations  in  F[qy,  ky) 
about  (ly  =  isin©  .  Therefore,  we  make  no  small-angle 
approximation,  either  in  our  far-field  analysis,  based  on 
(1 1),  or  in  our  solution  of  the  Schrbdinger  equation  for  the 
scattering  states  (3). 

3.  Device  Characteristics 

As  discussed  in  section  1,  ballistic  paths  on  the  order  of 
3  microns  have  been  obtained  in  GaAs/AlGaAs  HEMTs  at 
low  temperatures.  This  is  considerably  longer  than  the 
wavelength  of  a  I  meV  electron,  and  also  longer  than  the 
minimum  feature  sizes  achievable  by  electron  beam 
lithography  [16].  Because  of  this,  it  is  possible  to  choose 


the  dimensions  of  the  drain  region  in  Fig.  1  so  as  to  fulfill 
two  complementary  conditions:  on  the  one  hand,  the 
collectors  are  close  enough  that  a  large  fraction  of  the 
electrons  travel  ballistically  from  slit  gate  to  drain.  On  the 
other  hand,  they  are  far  enough  away  that  a  Fraunhofer 
diffiaction  pattern  is  detected.  This  can  be  achieved  by 
placing  the  tips  of  the  collectors  one  micron  from  the  slit. 
With  a  100  nm  pitch,  collectors  at  this  distance  would  yield 
an  angular  resolution  of  about  6  degrees. 

In  order  to  observe  diffraction  effects,  the  slit  width  must 
be  comparable  to  the  electron  de  Broglie  wavelength.  This 
in  turn  requires  a  treatment  of  the  states  that  is  based  directly 
on  the  Schrodinger  equation  (Section  2).  Having  placed  the 
collectors  in  the  Fraunhofer  regime  of  distances  tom  the  slit, 
however,  we  assure  that  the  potential,  which  depends  on  the 
source-drain  bias  Vds,  varies  with  a  characteristic  length 
scale  much  longer  than  the  wavelength.  This  allows  us  to 
retain  one  simplifying  assumption  of  the  usual  semiclassical 
approach  to  electronic  transport:  it  is  possible  to  neglect  the 
field  in  the  computation  of  the  eigenstates.  That  is,  we  can 
use  the  states  of  the  previous  section.  In  this  regime,  the 
dependence  of  the  current  on  Vos  arises  tom  variation  of  the 
state  occupation  probabilities. 

We  can  assign  an  occupation  probability  to  the  scattering 
states  (3)  by  observing  that  their  incident  current  density  is 
the  same  as  that  of  normalized  momentum  eigenstates  far 
from  the  slit  [11].  Furthermore,  we  observe  that  in 
equilibrium  (zero  bias)  there  is  no  net  current,  since  the 
current  transmitted  through  the  slit  from  scattering  states 
incident  from  the  source  at  any  energy  is  exactly 
compensated  by  a  deficit  of  reflected  current  associated  with 
degenerate  scattering  states  incident  tom  the  drain.  The  net 
current  density  in  the  drain  region  therefore  depends  only  on 
the  excess  occupation  probability  in  the  drain: 

/(0)  =jdky  jdk^  5n(k);(0;k) .  (12) 

0 

Hetey(0)d0  is  that  part  of  the  total  transmitted  current  which 
is  ejected  into  the  range  of  angles  (0,  0+d9),  The  current 
density  associated  with  incident  wave  vector  k, 

y(0;k)  =  — cos^0  k^  IF(tsin0,K)l^  ,  (13) 

Ttm 

is  computed  from  (1 1).  With  a  forward  bias  of  Vds,  the 
differential  occupation  probability  is 


6n(k)  =  n(£(k)-eVDs)-n(£(k))  ,  (14) 

with 

''(^)=l+exp(OT- 

We  have  evaluated  (12)  numerically.  Figure  2  shows  the 
result  fw  a  slit  width  of  250  nm,  an  electron  density  n  =  3  x 
10"  cm-t,  and  Vds  =  0.1  mV.  Lobes,  characteristic  of 
diffiaction,  can  be  seen  clearly. 

The  distribution  (14)  has  a  width  of  roughly  eVps  +  kT 
in  energy.  For  a  density  of  3  x  10"  cm-^  in  n-GaAs,  Ep  = 
11  meV,  compared  to  kT  =  0.36  meV  (at  4.2K).  The 
electron  distribution  can  therefore  be  made  approximately 
monochromatic  if  the  bias  is  kept  in  the  few-millivolt  range. 
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Fig.  2.  Angular  current  density  >(0)  for  a  slit  width  of  250 
nm  an  electron  density  n  =  3  x  10"  cm  ’,  and 
Vds  =  0.1  mV. 


Fig.  3.  Current  density  y(0)  for  a  slit  width  of  150  nm  and 
electron  density  n  =  3  x  10"  cm  ’,  for  biases 
Vds  of  (a)  3.16  mV.  (b)  1.0  mV,  (c)  0.316  mV  and 
(d)O.l  mV. 


and  for  V^s  in  this  range,  the  temperature  dependence  due 
to  the  width  of  (14)  is  negligible  (our  calculations  are 
performed  at  T  =  OK).  Thus,  obtaining  high  temporal 
coherence  is  not  difficult.  Curves  in  Fig.  3,  plotted  for 
different  values  of  the  b'as,  show  that  for  Vps  less  than 
about  1  mV,  tf ;  current  dcp  ity  is  essentially  proportional  to 
bias.  At  the  highest  bias,  the  diffraction  lobes  are  shifted 
toward  smaller  angles.  This  is  due  to  the  shorter  wavelength 
of  electrons  injected  at  high  bias.  (Calculations  in 
subsequent  figures  were  performed  using  Vos  =  0.1  mV.) 

The  above  suggests  that  there  is  only  a  weak  dependence 
on  temperature  below  about  40K.  In  the  real  device, 
scattering  will  broaden  any  sharp  features  that  (12)  may 
exhibit.  Nevenheless,  studies  using  hot  electron  injectors 
with  variable  injection  energy  [8]  display  a  pattern  that 
simplifies  our  analysis:  over  short  distances,  the  scattered 
(non-ballisdc)  component  of  the  electron  distribution  quickly 
assumes  a  broad  distribution  of  energies,  initially  very  flat 
over  the  whole  range  of  ener^es  from  zero  up  to  the  energy 
of  the  ballistic  electrons.  This  kind  of  scattered  component 
gives  rise  to  a  smooth,  featureless  background,  above  which 
die  diffraction  pattern  due  to  ballistic  electrons  should  be 
observable. 
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Another  source  of  incoherence  is  the  spatial  extent  of  the 
source.  A  common  measure  of  coherence  (appropriate  for 
quasi-monochromatic  sources)  is  Michelson's  fringe 
visibility  V  [10].  This  is  defined  in  terms  of  the  intensities 
of  an  adjacent  pair  of  extrema  in  the  diffraction  pattern  as 

V  =  /piK~ypuil  nsi 

JrtMxNmin 

The  most  commonly  studied  example  of  spatial 
coherence  effects  considers  a  double  slit  experiment,  where 
the  canonical  result  is  ihat  V  depends  on  a  quantity  ^afk, 
where  is  the  angle  subtended  by  the  source  as  viewed  from 
the  diffracting  slit,  and  equals  n  for  the  device  shown  in  Fig. 
1.  Diffraction  fringes  are  clearly  discernible  only  when  ([la^ 
is  small:  V  =  lsinc(7t^aA)l  (sinc(x)  s  sin(z)/x).  Applying 
this  classical  result  to  our  situation,  for  n  =  3x10"  cm  '  (Xp 
=  46  nm)  and  a  =  250  nm,  we  find  V  <  0.02.  In  fact,  this 
criterion  is  rather  conservative,  because  it  is  derived  in  the 
paraxial  approximation  of  small  6  and  relatively  large  aA. 
This  allows  one  to  ignore  the  disclination  factor  in 
Kirchoffs  rigorous  formulation  of  the  Huygens-Fresnel 
theory,  and  to  neglect  the  reduced  apparent  size  of  the  slit  as 
viewed  from  an  oblique  source.  The  effect  of  these 
neglected  factors  is  to  reduce  the  contribution  of  electron 
waves  incident  at  wide  angles,  and  effectively  focus  the 
incident  waves  into  a  source  distribution  with  smaller 
effective  extent  (ft. 

For  the  present  small  aA  situation,  y'm.x  varies  rapidly 
from  one  peak  to  the  next,  so  V  is  ambiguously  defined. 
For  a  conservative  estimate,  we  compare  the  current  density 
of  each  minimum  with  that  at  the  subsequent  maximum  (i.e., 
the  lower  of  the  two  adjacent  peaks).  In  Fig.  4  we  plot  the 
resulting  fringe  visibilities  for  the  lobes  of  Fig.  2,  and  find  a 
fairly  high  value  of  V  -  0.5.  It  is  unclear  at  present  to  what 
extent  these  higher  values  of  V  are  due  to  our  basic 
approximation  (6).  However,  it  seems  likely  that  the 
classical  value  0.02  may  be  regarded  as  a  lower  bound. 

4.  Device  Applications 

The  QUADFET  has  interesting  device  applications.  Two 
of  these  are  based  on  the  way  the  diffraction  pattern  varies 
with  gate  bias.  As  gate  bias  is  increased,  and  slit  width 
decreases,  the  diffraction  pattern  (Fig.  5)  and  its  lobes  (Fig. 
6)  shift  outward  to  larger  angles.  Viewed  at  a  particular 
collector,  this  leads  to  oscillations  in  the  current,  as  different 
lobes  cross  the  collector's  angular  position.  This 


Fig.  4.  Fringe  visibilities  V  for  the  250  nm  slit. 
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Fig.  5.  Current  densities  for  slit  widths  of  (a)  250  nm,  (b) 
200  nm,  (c)  150  nm  and  (d)  100  nm,  for  the  density 
and  bias  conditions  of  Fig.  1 . 


Fig.  6.  Positions  of  the  first  and  second  subsidiary  current 
maxima,  as  functions  of  the  slit  width. 


Slit  Width  (nm) 

Fig.  7.  Current  density  at  45'  viewed  as  a  function  of  slit 
width. 


phenomenon,  shown  in  Fig.  7,  is  essentially  an  oscillatory 
transconductance.  The  oscillatory  characteristic  makes 
possible  a  versatile  frequency  multiplier:  for  a  large  angle  or 
voluge  amplitude,  one  cycle  of  gate  voltage  can  correspond 
to  many  cycles  of  output  current  The  output  has  a  dominant 
component  at  an  adjustable  frequency,  unlike  ordinary 
ftequency  multipliers,  in  which  the  output  power  is  typically 
shared  by  a  range  of  harmonics  and  the  output  signal  must 


be  selected  by  filtering.  Since  the  device  is  ballistic  and  can 
operate  at  high  fiequencies,  muldplied  frequencies  can  range 
into  the  terahertz  under  proper  conditions  of  bias,  density 
and  angle. 

In  Edition,  collectOTS  can  be  chosen  which  are  turned  on 
or  off  in  complementary  response  to  the  input  voltage, 
allowing  the  formation  of  a  low-power  gate. 
Complementary  logic  may  thus  be  implemented  with 
QUADFETs  in  place  of  the  MOSFETs  used  in  CMOS. 
Possible  advantages  are  a  reduced  number  of  fabrication 
steps  and  the  inherent  high  speed  of  ballistic  transport. 

Furthermore,  QUADFETs  may  be  used  as  magnetic  field 
sensors.  A  magnetic  field  applied  perpendicular  to  the  plane 
of  the  device  curves  the  electron  trajectories,  shifting  and 
deforming  the  current  pattern  detected  at  the  drain  collectors. 


5.  Summary 

We  have  described  the  analysis  of  a  novel  field  effect 
transistor,  whose  operating  principle  is  the  quantum 
diffraction  of  electrons.  Appreciable  fringe  visibilities  can  be 
obtained  even  in  the  presence  of  low  source  coherence.  The 
useful  property  of  oscillatory  transconductance  was 
demonstrate  theoretically.  This  property  leads  to  a  number 
of  useful  device  applications. 
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